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Abstract 

We study a simple extension of the standard model to simultaneously explain neutrino masses, 
dark matter and the matter-antimatter asymmetry of the Universe. In our model, the baryon 
asymmetry is achieved by the leptogenesis mechanism, while the decaying dark matter with the 
lifetime of 0(10^^ s) provides a natural solution to the electron and positron excesses in Fermi 
and PAMELA satellite experiments. In particular, we emphasize that our model is sensitive to 
the structure at the endpoint around 1 TeV of the Fermi data. In addition, some of new particles 
proposed in the model are within the reach at the near future colliders, such as the Large Hadron 
Collider. 
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The observed neutrino oscillations and matter-antimatter asymmetry as well as the ev- 
idence for dark matter (DM) l| clearly imply physics beyond the standard model (SM). 
Recently, PAMELA [2] and ATIC [3] cosmic-ray measurements show the positron/electron 
excesses above the calculated backgrounds for the energy of (9(100) GeV. These data are 
consistent with the measurements of the high energy electrons and positrons fluxes in the 



cosmic ray spectra by PPB-BETS 



l], HEAT [si, AMS ^\ and HESS [?, 8]. Very recently, 
a more precise data by the Fermi LAT collaboration 9j also indicates some enhancements 
in the electrons + positrons flux in the 100 — 1000 GeV energy range. However, the Fermi's 
result is in conflict with the large excess in flux around 500 GeV range by ATIC. Similar 
conclusion has also been given by HESS based on the low energy data 8|]. In this study, we 
will concentrate on the combined data of PAMELA and Fermi (PF) without fitting that of 
ATIC. 

To explain the PAMELA/ATIC data, there have been many possible mechanisms, such 
as DM decays 10, 11|, DM annihilations [Q] and astrophysical sources 13 1, while recent 
studies related to the PF data without the ATIC one can be found in Refs. jlJ . fisj ]. In 
this paper, we would like to explore the possibility of connecting the neutrino masses to the 
dark matter problem as well as the baryon asymmetry of the Universe (BAU). In particular, 
we would like to pay attentions to models which can be tested directly by the future high 



16| and Linear Colhder (LC) [17 1 



energy colliders, such as the Large Hadron Collider 

We introduce three new neutral leptons Ni [i = 1, 2) and with the masses of Mi and 
M, and two new doublet scalars ( and rj with zero VEVs and the masses of and M^, 
respectively, in the SM. These new particles have non-trivial transformation properties under 
the two discrete symmetries Z2 and Z2 as listed in Table [H whereas the corresponding SM 
particles are trivial. 

TABLE L Transformations of the new particles under the discrete symmetries of Z2 and 
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Ni 


iV 


Z2 




+ 
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Z'2 


+ 




+ 





The relevant Majorana mass terms and Yukawa couplings as well as the soft breaking 



2 



term involving the new particles can be written as 

^^NfCN, + ^N^CN + y.MNj + y'MvN + liWC + H.c, (1) 

where i and j are the flavor indexes and Lj are the lepton doublets in the SM. We note 
that the soft breaking term in Eq. (1) breaks the two discrete symmetries to a diagonal 
one. In our study, we will assume the mass hierarchies of M^o < < Mi < M2 and 
< M < M^. We will demonstrate that the leptogenesis is achieved by A^^i decays, while 
is the decaying dark matter. 

The neutrino masses are generated by the one-loop diagram in Fig. [1] as proposed in 
Ref. 18| due to the quartic scalar interaction of 

^(0tC)2 + H.c., (2) 

where is the SM Higgs boson. A simple formula for the neutrino masses is found to 
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FIG. 1: Neutrino masses generated at one-loop level. 

be [18|, \li 
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(mujii = T / —V (3) 

k=l " 

with the SM Higgs VEV of v ^ 174 GeV. For the parameter set of A = C(10-^), = 
0(10-3) and Mi = 0(100 GeV - 10 TeV), we obtain = 0(0.01 - 0.1 eV), consistent 
with the current neutrino data. Similarly to the minimal seesaw model with two right- 
handed neutrinos 20|] , our model contains one massless neutrino with only normal or inverted 
hierarchy of the neutrino masses. However, the extended model with three Ni could still has 
the possibility of the quasi-degenerate neutrino masses. 



The leptogenesis mechanism in our paper is the same as that in Ref. 19j , provided by A^, 



and (. The relevant diagrams for the leptogenesis are shown in Fig. [2l The decay width of 




No 



FIG. 2: Tree level and one-loop diagrams for A'^i £i<^. 

Ni is given by 

mi - e^C^) = ^^M,r' (4) 
ibvr 



with r = 1 — M^/Mf. The out-of-equilibrium condition requires r ~ 10"'* 19|. We remark 
that this small value of r implies some degeneracy between Mi and M^. However, it may be 
avoided by including another new doublet with a soft breaking term similar to the discussion 
on the decaying dark matter discussed later. 

In the case of 3Mi < M2, the CP violating parameter in the leptogenesis is given by 

' ' -Im[(,t^)y^. (5) 



By using yij = O{10~^) and Mj = (9(100 GeV — 1 TeV) and assuming a maximal CP phase, 
the net BAU can be obtained as 

'^'B ^ ^ in-10 /c^ 

— ~ ~ 10 , (6) 

s 15 ' ^ ^ 

where g^, ~ 100 is the relativistic degrees of freedom. The decays of —>■ (^i'^u and 

help to avoid the dangerous relics from the singly- charged component of (, 

while C° may provide only a sub-dominant component of the DM due to the annihilation 



into gauge bosons 



2l|. 



The diagram for the DM decay is shown in Fig. [3l The decay width of is given by 



An \Mr,J M 



(7) 



where we have used the definition of 



2M 

By neglecting the effects of the off-shell C^, the lifetime of N is 



M± = (8) 



1 ttA'M , . 



FIG. 3: Diagram for the DM decay. 



with 



A 



(10) 



where we have included both charged and neutral modes. By taking y[ = 0(10^'*), yU = 
C(lkeV) and = C(100 TeV), one gets A ~ 10^^ leading to tn ~ lO^^ s with ~ 
0.5 TeV and M ~ 2 TeV. 

The normalized energy spectrum of the electron/positron for e^C^ can be written 

as 



dE 



2S{M_-E), 



while that for the decaying chain of — > C'^fi^{—>- e^2v) is 14 1 



dN, 



fie 



4 



(x3 - 1) - ^ (x^ - 1) 



subprocesses via exchanges of bosons, we have 



dN, 



dE 



dE 



dN^e 
dE 



dN^e 
dE '' 



where 



E„ 



dN, 



2Ml 



dE 



[m 



2 \2 



(m 



2 \2 



(11) 



(12) 



dE 3M_ ^ 

with X = E/M^ and < E < M_. For the electron/positron produced in the (^^ e^vQ^ 



(13) 



(14) 
(15) 
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with 



m,. 



El 



^ = {El + ElY-{\El\T\E, 



*|\2 



2m 



e; 



Ml 



m 



2m 



m^,^o = M^^ (1 



2E/M). 



(16) 
(17) 
(18) 



Note that we wiU concentrate on the case in which AM^ = — M^o is small to forbid 
the hadronic decay modes of W bosons. The normalized resultant energy spectrum of the 
electron/positron from the DM decays can be written as 



dN 
dE 



2 + e 



dN. 



dN 



+ ^~7Tr + ~1T^^ {Emax — E) 



(19) 



dE dE dE 

where e = |l/|jP/|l/eP have assumed 100% rate for the electron channel of the ( decay. 

We remark that we have neglected the tau-lepton effect in this study, but it can be included 
straightforwardly. 



The DM component of the primary electron/positron flux is given by 

M- 



10 



22 



dE'G{E, E'\ 



dN 

1e'' 



(20) 



where E is in units of GeV and c is the speed of light. All the information about astrophysics 
is encoded in the Green function of G{E, E'), given by 

10^6 



G{E,E') 



E^ 



exp[a + b{E'-^ - E"-'')]0{E' - E) [cm-^s] 



(21) 



where the normalization is adjusted to yield a local halo density p© ~ 1 GeVcm"^ 23|. 
We use the coefficients of a = —1.0203 and b = —1.4493 [22] for the spherically symmetric 
Navarro, Frenk and White density profile of the DM in our Galaxy 24] and the diffusion 



parameter 6 = 0.70 for the MED propagation model 25|, which is consistent with the 



observed Boron-to-Carbon ratio 



261 ]. Here, we have not taken account of the charge-sign 



dependent solar modulation [27| as well as other astrophysical uncertainties [25| , which could 
be significant in the energies below 10 GeV. The total electron and positron fluxes are 



K 



(22) 



6 



respectively, where is a primary astrophysical component, presumably originated from 



supernova remnants, $ _[+) is an exotic primary component from the DM decays, 



-(+) 



IS 



a secondary component from the spallation of cosmic rays on the interstellar medium, and k 
is a free parameter about 1 to fit the data when there is no DM primary source. We choose 



K = 0.7 to insure the flux calculation to be consistent with t 
fluxes, we will use the parameterizations obtained in Refs. 



27 



le data. For the background 
28| . given by 



11^0-9 + 3.2^2-15 

0.7^0-^ 



110^1-5 + 600^2-9 + 580^4-2 



[GeV^^cm-^s^^sr-^], 

[GeV~^cm-2s~^sr-^] 



1 + 650^2-3 + 1500^4-2 
where E is in units of GeV. 



[GeV-^cm-^s-^sr-i] 



(23) 
(24) 
(25) 
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FIG. 4: Electron + positron energy spectrum (left) and positron fraction (right) of the DM decays 
with TN = 2.5 X 1026 s, M = 2 TeV, = 500 GeV and AM^ = 1 GeV, where e = 1 (5) is 
represented by dashed (solid) lines, black points and blue rectangles stand for the observations of 
Fermi and HESS (left) and PAMELA and HEAT (right), and green triangles and dot-dashed lines 
correspond to the ATIC and backgrounds, respectively. 



In Fig. m [5], we show the electron plus positron energy spectrum (left) and the positron 
fraction (right) of the DM decays for = 2.5 x 10^6 [2 x lO^^] s, M = 2 TeV, = 500 GeV, 
AM^ = 1 GeV and e = 1 and 5, respectively, where the backgrounds are represented by dot- 
dashed lines. The electron + positron flux is multiplied by to compensate the roughly 
E'^ falling of the flux. From the figures, we see that the model with the enhanced muon 
effects is in good agreement with the Fermi, HESS, PAMELA and HEAT data. In particular, 
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FIG. 5: Legend is the same as Fig. 4 but with tn = 2 x 10^^ s. 

the energy spectrums in Figs. H] (left) and [5] (left) with e = 5 perfectly matches the Fermi's 
result. It is worth to mention that the results for the energy spectrum and positron fraction 
are not significantly dependent on the ( mass in the wide range of 100 GeV — 1 TeV besides 
the end point moving to a higher energy for a smaller M^. We remark that for a lighter DM 
particle, the drop in the electron flux occurs at a lower energy compared to the Fermi data. 



As for the collider signatures from the new partic 



es in the model, there are possible pair 



productions of ( directly by the SM gauge bosons 181]. However, the pair productions of 
Ni by the e~^e~ annihilations through the ( exchanges jisl are hard to be observed at the 
near future LC due to the small values of yij = O{10^^). In addition, although and 
1] will escape the detection in the next generation of colliders due to their weak couplings 
and heavy masses, their properties can be tested by precise measurements of the electron 
spectrum and positron fraction since the signals are not sensitive to the propagation models 
at the energies higher than 400 GeV. In particular, future measurements of the positron 
fraction at energies higher than 100 GeV can be crucial in testing the same origin of the 
Fermi and PAMELA electron and positron excesses. 

Finally, we remark that the neutral lepton A^ could be produced copiously after the big 
bang to become a typical unwanted relic. However, the inflation dilutes A^ away since its 
number density reduces exponentially. The present abundance of A^ may be generated by 
the e~^e~ annihilation. 

In conclusion, we have investigated a relatively simple extension of the SM to generate 
the small neutrino masses at one-loop level and the observed BAU by the leptogenesis 
mechanism. Our model also contains the decaying dark matter with the lifetime of O{10^^ s), 
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which provides a natural solution to the electron and positron excesses in cosmic rays in the 
energy range of 100 — 1000 GeV by Fermi and PAMELA. It should be emphasized that the 
structure at the endpoint around 1000 GeV of the Fermi data is crucial to determine the 
muon effects in the dark matter decays. More precise cosmic ray measurements around this 
energy range are clearly needed. 
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